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Thermodynamics and the phase structure of tire Polyakov loop-extended two flavor chiral quark- 
meson model (PQM) are explored. The analysis of the PQM model is based on the functional 
renormalization group (FRG) method. An appropriate truncation of the effective action with quarks 
coupled to background gluonic flelds is introduced. Within this scheme, we derive the renormal- 
ization group flow equation for the scale-dependent thermodynamic potential at flnite temperature 
and density in the presence of a symmetry breaking external field. The infiuence of fluctuations and 
of the background gluon field on the properties of net-quark number density fluctuations and their 
higher moments is explored. We study the dependence of the kurtosis of quark number fluctuations 
on the pion mass and show that, in the presence of a symmetry breaking term, the fluctuations lead 
to a smoothing of observables near the crossover transition. 

PACS numbers: 



I. INTRODUCTION 

The phase diagram of strongly interacting matter at 
nonzero baryon density and high temperature has been 
a subject of growing interest in recent years. Calcu- 
lations done within the Lattice Gauge Theory (LGT) 
show a clear separation between the confincd-hadronic 
and deconfined, quark-gluon plasma phase at finite tem- 
perature Quantum Chromodynamics (QCD) ex- 
hibits both dynamical chiral symmetry breaking and con- 
finement at finite temperature and densities. However, 
since QCD thermodynamics at large baryon densities is 
still not accessible for first principle LGT calculations, 
many phenomenological models and effective theories 
have been developed. 

The hadronic properties at low energy as well as the 
nature of the chiral phase transition at finite tempera- 
ture and densities have been successfully described and 
explored in such effective models. Also, the physics of 
color deconfinement and its relation to chiral symmetry 
braking has been recently studied in terms of effective 
chiral models [2l4l3j . The idea to extend the existing chi- 
ral models, such as the Nambu-Jona-Lasinio (NJL) or 
the chiral quark-meson model, by introducing couplings 
of quarks to uniform temporal background gauge fields 
(Polyakov loops) was an important step forward in these 
studies 

It was shown that the Polyakov loop extended Nambu- 
Jona-Lasinio (PNJL) or quark meson (PQM)[ll| 
models reproduce essential properties of QCD thermody- 
namics obtained in the LGT calculations already in the 
mean-field approximation. However, to correctly account 
for the critical behavior and scaling properties near the 
chiral phase transitions one needs to go beyond the mean- 
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field approximation and include quantum fluctuations 
and non-perturbative dynamics. This can be accounted 
for by using the functional renormalization group (FRG) 
UJ-UJ. Until now this method was applied in the NJL 
and quark-meson model, where the FRG equation was 
formulated for quarks coupled to meson fields jl8l - l2^ . 

In this paper we propose a truncation of the PQM 
model which is suitable for a functional renormaliza- 
tion group analysis. We introduce an additional cou- 
pling of the chiral condensate to a gluonic background via 
Polyakov loops. In this way the Polyakov loop dynam- 
ics is represented by a corresponding background tempo- 
ral gauge field. We thus use the functional renormaliza- 
tion group approach in order to include fluctuations of 
the meson fields, while the Polyakov loop is treated on 
a mean-field level. Consequently, our calculation lacks 
an FRG description of the fluctuations of the Polyakov 
loops. However, already in this approximation we flnd 
an important effect of the interactions of quarks with 
the effective gluon fields on the thermodynamics. For 
comparison we also present results for the chiral quark- 
meson model (QM) without Polyakov loops in the FRG 
approach in order to make the difference in the thermody- 
namics of the two effective models more transparent and 
to emphasize the role of the gluonic background field. We 
also compare with the PQM and QM models in mean- 
field approximation. 

In this study we use the Taylor expansion method to 
describe the thermodynamics at finite chemical poten- 
tial /i. We discuss the influence of fluctuations and the 
background gluon field on the properties of the net-quark 
number density susceptibilities. Results on the proper- 
ties of generalized susceptibilities of the net quark num- 
ber density in the presence of mesonic fluctuations in a 
gluonic background are presented. We also explore ways 
to reduce the dependence on the ultraviolet cutoff in the 
FRG approach. 

The paper is organized as follows: In Sec. II we in- 
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troduce the PQM model. In Sec. Ill we formulate the 
renormalization group equation for the thermodynamic 
potential in the PQM model. In Sec. IV we consider the 
mean-field approximation. In Sec. V we discuss the crit- 
ical properties and thermodynamics of the PQM model 
in the presence of mesonic fluctuations. 



II. THE POLYAKOV-QUARK-MESON MODEL 

The chiral quark-meson model is used as an effective 
realization of the low-energy sector of the QCD. How- 
ever, because the local SU (Nc) invariance of QCD is re- 
placed by a global symmetry, one cannot describe de- 
confinement phenomena in this model. Recently, it was 
argued that, by connecting the chiral quark-meson model 
with the Polyakov loop potential, the confining properties 
of QCD can be approximately accounted for [HI. 
Consequently, the Polyakov-quark-meson (PQM) model 
effectively combines both the chiral and confining prop- 
erties of QCD. 

The Lagrangian of the PQM model reads 
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C = q[ip- g{(J + h5T7r)] q + ^{d^af + -{d^nf 



-[/(a,7f) -U{£J*) 



(1) 



The coupling between the effective gluon field and the 
quarks is implemented through the covariant derivative 



da - iA. 



(2) 



^ — The spatial components of the 



where A 

gluon field are neglected, i.e. A^ ^ ^^o^o- Moreover 
U {£,£*) is the effective potential for the gluon field ex- 
pressed in terms of the thermal expectation values of the 
color trace of the Polyakov loop and its conjugate 



N, 



■(Trei(x)), r 



iTrcLHx)), (3) 



with 



L{x) = V exp 



I I dTAi{x^T) 




(4) 



where V stands for path ordering, /3 = I/T and = 
IAq. 

The 0(4) representation of the meson fields \s cf) ~ 
{a, tt) and the corresponding SU (2) l®SU (2)^^ chiral rep- 
resentation is given hy a + it ■ 7175 . There are Nf = 4 
mesonic degrees of freedom coupled to Nf = 2 flavors of 
quarks. 

The purely mesonic potential of the model, U (a, tt ) , is 
defined as 



C/(a,7f) 



A 



ccr. 



(5) 



while the effective potential of the gluon field is 
parametrized to preserve the Z{5) invariance: 

^^i^ = J^te-^-^if + t^) + ^-^{t£)\ (6) 

The parameters, 

hiT) = ao+ar(f)+a,(^)' + a3(f)' (7) 

with ao = 6.75, ai = -1.95, as = 2.625, 03 = -7.44, 
63 = 0.75 and 64 = 7.5 were chosen to reproduce the 
equation of state of pure gauge degrees of freedom calcu- 
lated on the lattice. At the temperature Tq = 270 MeV, 
the critical temperature obtained for pure gauge theory, 
the potential ^ yields a first order phase transition. 



III. THE FRG METHOD IN THE PQM MODEL 

The functional renormalization group is an important 
tool for addressing nonpcrturbative problems within the 
quantum field theory. It is based on an infrared (IR) reg- 
ularization with the momentum scale parameter, k, of 
the full propagator which turns the corresponding effec- 
tive action into a scale dependent functional [Tj-[l3]. 
The change of Pfe with the change of the momentum scale 
is described through the flow equation: 



where r[,^'°'' and P^^'^'' denote the second fimctional 
derivative of Pfc[$,'0] with respect to the bosonic ($) 
and fermionic (■0) fields, respectively. These derivatives 
correspond to the inverse of the full propagators at the 
scale k. The trace in Eq. ([5]) denotes a momentum in- 
tegration and a summation over all internal indices (e.g. 
flavor, color, and/or Dirac). Here $ and ip denote bosonic 
and fermionic fields, respectively. The effective average 
action, P^, governs the dynamics of a theory at a mo- 
mentum scale k and interpolates between the bare ac- 
tion, P/c^A = S, and the full quantum effective action, 
Tk=Q = P. The regulator function, Rk, describes how the 
small momentum modes are cut off and is to some extent 
arbitrary (l7j . The derivative of this function, dkRk, im- 
plements the Wilsonian idea of successively integrating 
out momentum shells. 

In the PQM model the formulation of the FRG flow 
equation ^ would required an implementation of the 
Polyakov loop as a dynamical field. However, in this 
exploratory calculation we treat the Polyakov loop on 
the mean-field level. This allows us to formulate the FRG 
flow equation for the truncated effective PQM action in 



Tr {dkRkF (Pi,°'''[$,V] + i?fcF 



.(0,2) 



(8) 
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Euclidean space-time {t 



potential, Aq, are introduced through the following sub- 
stitution of the time derivative 



+ gq{a + iT- TT-fr,) q + Uk{p)} , 
where the field p is given by 



(9) 



dr ^ dr — (/i + iAi). 



(11) 



Since the gauge potential, Aq, is treated as an effec- 
tive background gluonic field, which does not flow with 
the scale parameter fc, the explicit contribution of the 
(10) Polyakov loop potential to the truncated effective action 
is suppressed at this stage. Later on we will restore the 
The finite quark chemical potential, fi, and the gauge Polyakov loop potential. 
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As in the previous studies j25| . we employ the optimized regulator functions [26| in the numerical implementation. 
For bosons, this cutoff function depends only on the spatial components of the momentum. 



whereas for fermions we use [25[ 



J [qo + la^Y + q2 j 



(12) 



(13) 



where = p + iAo is modified due to the coupling to the background gluon field. 

Using Eq. ([9|), with the cutoff functions given above, and the relation Q = TT (2^, we obtain the flow equation 
for the scale dependent grand canonical potential for the quark and mesonic subsystems 
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l + 2nBiE,;T) 



(14) 



l~N{£,t;T,p)-Ni£,£*;T,p) 



Here nBiE^^^^'jT) is the bosonic distribution function 

nB{Eyr,a;T) = 

with the pion and sigma energy 

E^ = \/k^ + nl , E, 



1 



eME.,JT) - 1 



/c2 + fifc + 2pSlfe 



where the primes denote derivatives with respect to p and = fi + ccr. The functions N{£, £*; T, p) and N{£, £*;T, p), 
defined by 



...^ . ^ 1 + 21* exp[/?(£;, - p)] +£expmE, - p)] 

^ ' ' l + M exp[2p{Eg ~ p)] + U* exp[/3(£;g - p)] + exp[3;3(£;5 - p)] ' 

N{£,t-T,p) = N{t\£-T,~p), 

are fermionic distributions, modified due to the coupling to gluons, and 



E, = ^k^ + 2g\ 



P 



(15) 
(16) 

(17) 



is the quark energy. In the absence of the background gluon field, e.g. in the QM model, £,£* 1 and the usual 
Fermi-Dirac distribution function is recovered; N(\,1]T, p) — np(E„;T, p) = l/{exp[(£'„ — p)/T] + 1}. 
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In order to solve the flow equation (jl4p . we expand Clk{T,p)^ in a Taylor series around the scale-dependent 
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minimum at cr^ 



(18) 



Within the FRG approach, the thermodynamic potential 
of the QM model is obtained from Eq. ([20|) . by dropping 
the effective Polyakov loop potential and setting I ~ i* ~ 
1. 



As in Ref. [25|, we truncate the expansion at m = 3. 
The Taylor coefficients, Oi^fe, are functions of the scale, fc, 
the temperature, T, the chemical potential, fj,, and the 
Polyakov loop, £ and £* . 

The minimum of the thermodynamic potential is de- 
termined by the stationarity condition 



da 



da 



-c = 0. 



(19) 



The set of flow equations for each expansion coefficient, 
Ui^ki arising from the Eq. (jlSp is solved numerically with 
an initial cutoff A = 1.2 GeV. The initial conditions, 
Q.i.fc=A and c, are chosen to reproduce the physical pion 
mass TOtt = 138 MeV, the pion decay constant /^r = 93 
MeV, the sigma mass mo- = 700 MeV and the constituent 
quark mass mg = 335 MeV at the scale fc = in vacuum. 
The symmetry breaking term, c = rn^f.,^, corresponds to 
an external field and consequently does not flow. In this 
work, we neglect the flow of the Yukawa coupling, g. 

By solving equation (fTS]) we obtain a thermody- 
namic potential for the quark and mesonic subsystems, 
^k^o{(-: i*;T, /x), as a function of Polyakov loop variables 
£ and £* . So far these variables are arbitrary. The full 
thermodynamic potential in the PQM model, including 
quark, meson and gluon degrees of freedom is obtained by 
adding to nk^o{£,£*;T,fi) the effective gluon potential 
U {£,£*) from Eq. ©, 



ni£, r ; T, ^l) = n^^oii, t-r, ^) + u{£, t). 



(20) 



At a given T and the Polyakov loop variables, £ and 
I* are determined by the stationarity conditions: 



^r!(£,r;T,/.) = 



0, 



(21) 
(22) 



J 



In the following we explore the thermodynamics of the 
PQM-FRG model. The role of the gluonic sector is as- 
sessed by juxtaposing the results obtained within the 
PQM and QM models, while the effect of mesonic fluc- 
tuations in the PQM model are deduced by comparing 
with the mean-field approximation. 



IV. THE MEAN-FIELD APPROXIMATION 



The importance of mesonic fluctuations is best illus- 
trated by comparing the FRG approach with the mean- 
field (MF) approximation for mesons. In the latter, both 
quantum and thermal fiuctuations are neglected and the 
mesonic fields are replaced by their classical expecta- 
tion values. In the formulation presented in the previ- 
ous section, the mean-field approximation is recovered 
by omitting the pion- and sigma-meson contributions in 
Eq. (|T4l) . Alternatively, this approximation can be ob- 
tained directly from the partition function for the La- 
grangian ([T]) by integrating out the quark degrees of 
freedom. These approaches are equivalent up to a sur- 
face term, which is negligible for large cuttoffs, A — > oo. 
We explicitly retain the divergent vacuum contribution, 
which is regularized by an ultraviolet cuttoff Aj\/i? . 
The importance of this contribution was demonstrated 
in Refs. fHIH. 



A detailed derivation of the mean-field approximation of PMQ model can be found in Ref. [3C 
therein. Here we present only the final result for thermodynamic potential 

^Imf =U{l,t) + U{{a),{n) = Q) + nqq{{a), £,£*). 

The contribution of quarks with mass = g{a) is given by 



and references 



(23) 



n,^{{a),£,£*) = -2NfT 



d^p {N^Eg^^^^ 



(2^)3 



e{Al,p -p^)+\n g^+\{a) T, /.) + In g^-\{a) , T ; T, /i) , (24) 



where 



gW((a),£,r;r,/i) = I + Si e^-iEg - ij)/T] 
g(-)((a),Ar;r,/i) = g^+H{a) , £* J;T, -fi); 



M* exp[-2{Eg - ^i)/T] + exp[-3(£;, - /i)/r], 



(25) 
(26) 
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FIG. 1: Thermal average of the Polyakov loop, i, and of the order parameter of the chiral phase transition, (a), as functions of 
temperature at zero baryon chemical potential in the mean- field approximation(left panel) and in FRG approach (right panel). 
The solid and dashed lines are obtained in the PQM model, while the dash-dotted lines correspond to the QM model. 
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FIG. 2: Temperature derivatives of thermal average of the Polyakov loop, and of the order parameter of the chiral phase 
transition, (a), as functions of temperature at zero baryon chemical potential in the mean-field approximation(left panel) and 
in the FRG approach (right panel) . The notation is the same as in Fig. [1] 



and Eq — 



4- is the quark quasi-particle energy. 



The equations of motion for the mean fields are obtained 
by requiring that the thermodynamic potential be sta- 
tionary with respect to changes of cr, I and t*: 



MF 



MF 



MF 



da 



dl 



dt 



(27) 



The model parameters are fixed to reproduce the same 
vacuum physics as in the FRG calculation, as described 
in the previous section. The additional free parameter, 
AmFi is chosen so as to reproduce the quark condensate 
in vacuum (uu) = —(260 MeV)^; we find Amf — 674 
MeV. 



V. FLUCTUATIONS AND 
THERMODYNAMICS OF THE PQM MODEL 



The thermodynamic potential obtained in section IIIII 
can be used to explore the influence of the gluonic back- 
ground field on the thermodynamics in the PQM model 
including the effect of fluctuations. In previous studies of 
the chiral quark-meson model within the FRG approach, 
it was shown that non-perturbative meson contributions 
modify the position of the chiral boundary and the criti- 
cal end point (CEP) in the (T, /x)-plane. The pseudocrit- 
ical temperature and chemical potential is usually identi- 
flcd by a maximum in the temperature derivative of the 
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order parameter or in the chiral susceptibihty. In Figs. [T] 
and [5] we show the dependence of (cr) and Polyakov loop 
I and of their temperature derivatives computed in the 
FRG approach as well as in the mean-field approximation 
for the PQM and QM models. 

As expected, for a physical pion mass the model ex- 
hibits a smooth crossover as a function of temperature. 
Moreover, as shown in Fig. [2J the temperature deriva- 
tives of both the Polyakov loop and the chiral order 
parameter exhibit peak structures at approximately the 
same temperature, = 257 MeV in the PQM-FRG 
and Tc = 230 MeV in the PQM-MF model. In general 
the relation between the deconfinement and chiral phase 
transitions is not well established. Recent LGT calcula- 
tions lead to conflicting conclusions concerning the rela- 
tive position of these transitions [ll. |3]|. 

A comparison of the QM and PQM models shows that 
the inclusion of gluon degrees of freedom shifts the chi- 
ral phase transition to higher temperatures both with 
and without mesonic fluctuations. The latter lead to a 
smoothing of the transition, as shown in Figs. [1] and [2] 



A. Quark number fluctuations 

To explore the influence of non-zero baryon chemical 
potential on the thermodynamics we employ the Taylor- 
expansion in ^/T around /i = applied in LGT [l| and 
in model calculations [13, H^- In particular, we obtain 
the expansion of the thermodynamic pressure 



where 



Cn{T) 



OO -| 



n=0 



9"[p(r,/.)/r4 



^{^l/Ty 



(28) 



(29) 



The expansion coefficients Cn{T) are generalized sus- 
ceptibilities that characterize the fluctuations of the net 
quark number 5Nq = Nq—{Nq) at vanishing chemical po- 
tential [Ui36i. In particular, the first two non- vanishing 
derivatives, C2 and C4, are the second and fourth order 
cumulants: 



C2 = || = ((<5iV,)^), 
C4 - {{SNq)^) - 3{{SN, 



(30) 



where Xq is the regular quark number susceptibility. 

The temperature dependence of the coefficients C2 and 
C4, obtained in the mean- field approximation and in the 
FRG approach, is shown in Figs. [3] and 21 The coeffi- 
cient C2 increases monotonously with temperature. In 
the mean-field approximation, C2 increases rapidly in the 
critical region and approaches the ideal gas limit at high 
temperatures. The FRG results show a rather different 



behavior at high temperatures. The corresponding sus- 
ceptibility is strongly suppressed above the chiral tran- 
sition, due to the ultraviolet cutoff. The contribution 
of momenta beyond the cutoff to the thermodynamics is 
missing. In order to obtain the correct high-temperature 
behavior of the susceptibilities and other observables, one 
may supplement the FRG with the perturbative contri- 
bution of the high-momentum states. A simple model 
for implementing this correction was proposed in Ref. 
[S^l, where the corresponding contribution {k > A) to 
the fiow is approximated by that of a non-interacting gas 
of quarks and gluons. The flow equation for the high 
momentum contribution to the QM model then reads 



^^{2(iV2-l)[l + 2nB(fc;T)] 
1 - nF(fc;T, /i) - nF(fc;T, -/i) |, (31) 



127 



where the dynamical quark mass is neglected, i.e. we 
set Eq = k. For the PQM model we generalize this pro- 
cedure, by including the interaction of quarks with the 
Polyakov loop 



NcNfk^ 
37r2 



(32) 



l-Nie,t;T,ti)~Nie,t;T,fi) 



Since the effective gluon potential, U {£,£*), is fitted to 
reproduce the Stefan Boltzmann limit at high tempera- 
tures, the explicit gluon contribution is omitted for con- 
sistency. 

We thus proceed as follows: Eq. ((3T|) or Eq. ((32)) is 
integrated from fc = 00 to fc = A, where we switch to 
the QM or PQM flow equation, Eq. (HH). The divergent 
terms in the high-momentum flow equations pill32p are 
independent of mesonic and gluonic fields, and of temper- 
ature and chemical potential. Consequently, they can be 
absorbed in an unobservable constant shift of the ther- 
modynamic potential. 

In the right panels of Figs. |3] and |3]wc show the results 
for C2 and C4 in both the QM and PQM models. We find 
a smooth transition between the high and low energy 
regimes of the theory. The Stefan Boltzmann limit, C2 = 
2, is reproduced at high temperatures. 

In contrast to C2, the coefficient C4 shows a non- 
monotonic behavior near the chiral pseudocritical tem- 
perature. As shown in Fig. HI C4 exhibits a pronounced 
peak at the pseudocritical temperature. The shape of 
the peak is changed considerably when mesonic fluctua- 
tions arc included; its maximum value is reduced while 
the width is enhanced. We note that the peak value ob- 
tained in the mean-field approximation is strongly depen- 
dent on the cutoff Amf, used to regularize the divergent 
vacuum contribution. A smaller value of the cutoff yields 
a larger maximum value of C4. The significance of the 
fermion vacuum contribution was recently emphasized in 

Refs. (H,!!!. 
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FIG. 3: The coefficient C2 as a function of temperature at zero baryon cfiemical potential for the PQM and QM models in the 
mean-field approximation (left panel) and in the FRG approach (right panel). The thin lines in the rrightpanel indicate the 
FRG results after inclusion of the high-momentum flow (see text for details). 
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FIG. 4: The coefflcient C4 as a function of temperature at zero baryon chemical potential for the PQM and QM models in the 
mean-field approximation (left panel) and in the FRG approach (right panel). As in Fig. |3] the thin lines indicate the FRG 
results after inclusion of the high- momentum flow (see text for details). 




FIG. 5: The kurtosis 7?4,2 as a function of temperature at zero baryon chemical potential for the PQM and QM models in the 
mean- field approximation (left panel) and in the FRG approach (right panel). The thin lines indicate the FRG results after 
inclusion of the perturbative high-momentum modes (see text for details). 
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The general trends of the mean-field calculation can 
be inferred from Landau's theory of critical phenomena, 
where the thermodynamic potential is assumed to be a 
polynomial in the order parameter a 

n{T, ^i; a) = n^giT, ^i) + ia(r, ^i)a^ + ha" + ca. (33) 

Here flhg{T,fi) is a background contribution, indepen- 
dent of cr, and c is a symmetry breaking term. The co- 
efficient a is generally assumed to be a linear function 
of the temperature, a = A ■ (T — Tc), where Tc is the 
critical temperature of the second-order phase transition 
for c ~ 0. An extension to small, non-zero values of the 
chemical potential is achieved by assuming 

a{T,^i) = A-{T-T,) + BiJ.^, (34) 

where both coefficients A and B are positive. In general 
the effective quartic coupling constant, 6 > 0, also de- 
pends on temperature and chemical potential. However 
this dependence is irrelevant for T ~ Tc and ^ ~ 0. In 
the chiral limit, c = 0, the C2 and C4 coefficients are easily 
obtained from Eq. ([55]) : 

A ^ ^ 

C2 - c,,g + —^r^0{T,^T), (35) 
C4 = Cii,g + ^OiT,-T). (36) 

Thus, C2 is not differentiable at the critical point, while 
C4 exhibits a discontinuity. The background parts, C2bg 
and Cibg, are smooth functions of temperature and do 
not change the critical behavior. For a finite pion mass 
(c 7^ 0), the transition is of the cross over type and the 
sharp structures in C2 and C4 are smoothened. Conse- 
quently, in the QM and PQM models, the peak structure 
appearing in C4 is due to the decreasing quark mass, and 
thus closely connected to the restoration of chiral sym- 
metry. By contrast, in the resonance gas model, the peak 
in C4 is due to the increasing contribution of higher res- 
onances. 

In general, the singular part of thermodynamic poten- 
tial is in chiral limit controlled by the critical exponents 
of the tree-dimensional 0(4) symmetric spin model. Con- 
sequently, at zero chemical potential: 

ns^ng ~ {T - T,f~" , (37) 
C2s^ng ~ (T-Tj^-", (38) 
C4s^ng ^ (T-Tj-". (39) 

Critical fluctuations renormalize the exponent a from its 
mean-field value q;mf = to a = —0.21. Consequently, 
fluctuations lead to a weakening of the singularity in the 
chiral limit. This is reflected in a smoothening of the 
temperature dependence for finite pion mass, as seen in 
Figs. [3] and m We note that in QCD the singularity in C2 
and C4 may be obscured by resonances, which presumably 
dominate the thermodynamics of the low-temperature 
phase. 



A non- vanishing background gluon field leads to qual- 
itative changes of the quark fluctuations. In the low- 
temperature phase, there is a strong suppression of the 
quark fluctuations due to their coupling to the Polyakov 
loop (see Figs. [3]and|4|). This is because single and dou- 
ble quark modes are suppressed, when the expectation 
value of the Polyakov loop, i and i* , is small (see Eqs. 

(fniT5|) and mm ). 

Lattice studies of QCD at finite temperature and den- 
sity as well as chiral model calculations show, that the 
ratio (kurtosis) 

i?4 2 = - (40) 

C2 

is an useful probe of the deconfinement and chiral phase 
transitions [3^ [s^ [l^ . In the high and low temperature 
regimes, the kurtosis refiects the net quark content of 
the dominant baryon number carrying effective degrees 
of freedom [H, [s^ . 

In Fig. [5] we show the ratio _R4.2 as a function of tem- 
perature. Both in the FRG approach and in the mean- 
field approximation for the PQM model, the kurtosis 
drops from i?4^2 — 9 to i?4^2 < 1 in the transition re- 
gion, as expected due to the change in quark content of 
the baryon carrying effective degrees of freedom . At 
low temperatures the effective three quark states dom- 
inate, while at high temperatures single quarks prevail. 
In the mean-field approximation (see the left panel of 
Fig. [5]), the kurtosis exhibits a well defined peak at the 
transition temperature. The height of the peak depends 
on the value of the pion mass [H, . This dependence 
is illustrated in Fig. [5] both in the mean-field approxima- 
tion and in the FRG approach. The inclusion of mesonic 
fluctuations weakens the dependence on the pion mass. 
Thus, for a physical value of mjr, the kurtosis decreases 
monotonously with temperature in the FRG approach. 
A peak in i?4^2 appears only for a pion mass below its 
physical value. In the mean-field approximation, the de- 
pendence on the pion mass is more pronounced, but still 
not as dramatic as found in Ref. [33|. This is due to 
the different treatments of the divergent fermion vac- 
uum contribution. In Ref. [38j the vacuum term was 
subtracted at T = and its dependence on temperature 
was neglected. Here we regularize the vacuum term with 
an ultraviolet cut off and retain its temperature depen- 
dence, which is due to the in-medium quark mass. The 
role of the vacuum term in mean-field calculations and 
its influence on the chiral phase transition are explored 
in detail in Ref. 29]. 

In effective quark models where the Polyakov loop 
is neglected, the thermodynamics is governed by one- 
quark states at all temperatures. Consequently, in the 
QM model the low-temperature limit of the kurtosis is 
i?4^2 = 1 rather than 9. Consequently, in such models, 
the interaction of quarks with the background Polyakov 
loop is essential for reproducing the low-temperature be- 
havior of i?4_2 found in LGT calculations. 
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VI. SUMMARY AND CONCLUSIONS 

We have discussed thermodynamic properties of the 
Polyakov loop extended quark-meson effective chiral 
model (PQM), including fluctuations within the func- 
tional renormalization group method (FRG). A trunca- 
tion of the PQM model was introduced, which allowed us 
to extend previous renormalization group studies by in- 
troducing the coupling of fermions to the Polyakov loop. 
We have thus formulated and solved the flow equation for 
the scale dependent thermodynamic potential at finite 
temperature and density in the presence of a background 
gluonic field. 

In our studies of the thermodynamics we have dis- 
cussed the role of fluctuations on the critical properties 
of the model. The fluctuations of the net quark num- 
ber density and the ratio of the fourth to second order 
cumulants, the kurtosis, were analyzed. 

Furthermore, a comparison of the FRG approach with 
and without a gluon background field was performed. 
The influence of ultraviolet cutoff effects on the thermo- 



dynamics within the FRG approach was also discussed. 
We have shown that the FRG extended quark-meson 
model preserves the basic properties of the kurtosis, ob- 
tained in LGT calculations only if the Polyakov loop is 
included in the flow equation. Thus, the extension of the 
FRG method proposed here, accounting for the coupling 
of fermions to background gluon flelds, is of particular rel- 
evance for effective descriptions of QCD thermodynamics 
near the phase transition in terms of the quark-meson 
model. 
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